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Abstract. Using the Yguazu-a three-dimensional hydrodynamic code, we have computed a set of numerical simulations of 
heavy, supersonic, radiatively cooling jets including variabilities in both the ejection direction (precession) and the jet velocity 
(intermittence). In order to investigate the effects of jet rotation on the shape of the line profiles, we also introduce an initial 
toroidal rotation velocity profile, in agreement with some recent observational evidence found in jets from T Tauri stars which 
seems to support the presence of a rotation velocity pattern inside the jet beam, near the jet production region. Since the Yguazu- 
a code includes an atomic/ionic network, we are able to compute the emission coefficients for several emission lines, and we 
generate line profiles for the Ha, [O I]/16300, [S II]/16716 and [N II]/16548 lines. Using initial parameters that are suitable for 
the DG Tau microjet, we show that the computed radial velocity shift for the medium-velocity component of the line profile 
as a function of distance from the jet axis is strikingly similar for rotating and non-rotating jet models. These findings lead us 
to put forward some caveats on the interpretation of the observed radial velocity distribution from a few outflows from young 
stellar objects, and we claim that these data should not be directly used as a doubtless confirmation of the magnetocentrifugal 
wind acceleration models. 
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■ 1 . Introduction disk in which the magnetic field lines have their origin in the 

disk itself. They analyzed cold and warm models. These mod- 
It is now widely accepted that Herbig-Haro jets are a common els are mled by dimen sionless parameters: the disk aspect 
by-product of the formation of a low mass star. In itially, it was ratio s = hjr ^ where h is me disk scale height at me cylindrical 
proposed that these je ts originate as stellar winds (|Cant6 1980| radius r> the turbulence level a m = v m /V A h (where v m is the 
|HartmannetaL_I982). However, the apparent correlation be- turbu i ent ma g ne tic diffusivity and V A is the Alfven velocity at 
tween ejection and accretion signatures favoured disk wind the disk mid p lan e), and the ejection index £ = d In M a / d In 
models. Blandford & Payne f[982 1 proposed a model in which r (where ^ a is the accretion rate), which gives the efficiency 
a wind is launched from the disk surface through open mag- of the launching proce ss. In the cold models, the enthalpy at 
netic field lines which make an angle greater than 30° with the jet basis is nulL while the warm jet launching models have 
the vertical direction. This was called the magneto-centrifugal some heating source that provides additional thermal energy, 
model because the centrifugal force exerted in a particle at- xhe ej ection efficiency for cold jets is at most £ ~ 0.01 while 
tached to an open magnetic field line (satisfying the condi- war m jets could have efficiencies of up to £ ~ 0.5. Another kind 
tion mentioned above) accelerates it to speeds above the es- of disk model is the so called x _ wind (Shu et a l.|T994l In this 
cape velocity. More rec ently, Ferreira 11922) and Casse & case> the disk is exact i y at the corotation point and 
Fen-eira ( 2000a 2000b I proposed self-similar jet launching the jet is i aun ched from this truncation radius. In both models, 
mechanisms arising from an extended region of the accretion the rotation of the disk plays a very important ro i e . The exact 

point or region from which the jet arises is still under debate. 
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Although the magneto-centrifugal acceleration (hereafter, 
MCA) is widely accepted as being the general mechanism be- 
hind the launching of a disk wind, the lack of observational 
support for it until recently has contributed to keep the sub- 
ject open. In the last few years, however, observations of a few 
jets associated with young T Tauri stars appear to support this 
model, since the interpretation of the radial velocity profiles 
close to the launching region is both qualitatively and quan- 
titatively in agreement with them (e.g., Bacciotti et al. 2002 
Coffey et al. [2004 Woitas et al. 2005 1. These authors claim 
that, for a few jets associated with T Tauri stars (namely DG 
Tau, Th28, RW Aur and LkHa321), there is a trend in the 
radial velocity pattern which indicates that the jet is rotating. 
Furthermore, they find that the rotation velocity implied by the 
observed radial velocity distribution is in good agreement with 
the rotational profile expected for a wind launched by the MCA 
model (see Pesenti et al. 2004 1. 

Bacciotti et al. ( 2002 1 have carried out HST/STIS measure- 
ments of the [O I]^6300,6363 and [S II]X16716,6731 emis- 
sion line profiles of the DG Tau jet. The observation was de- 
signed to construct a datacube, changing the position of the 
long slit (of a/1 x 52'.') several times. The datacube was then 
used to determine the spatial distribution of the radial velocity. 
In particular, they have taken spectra for seven different posi- 
tions perpendicular to the DG Tau jet, for a given distance from 
the source, by offsetting each slit position by 0707. Assuming 
that the jet is axisymmetric, they carefully determine the po- 
sition of the jet axis, by determining the peak of each emis- 
sion line, and radial velocities associated with symmetrically 
disposed positions (on both sides of the jet axis) were mea- 
sured. The line profiles obtained by Bacciotti et al. (2002) were 
decomposed into two components, and the medium-velocity 
component (hereafter, MVC) was used to perform the analy- 
sis. The radial velocity (associated with the MVC) was then 
plotted against the slit position (i.e., distance from the jet axis), 
for four regions of increasing distance from the source (named 
regions I, II, III and IV, at distances from the source of 07075, 
07175, a/275 and a/4, respectively). They found that the radial 
velocity shift for positions symmetrically distant from the jet 
axis is systematically negative, regardless of the chosen region 
(I, II, III or IV) or the emission line 1 . This systematic negative 
shift in the radial velocity differences was interpreted by the 
authors as evidence that the jet is rotating. 

The values found for the rotational velocity of the DG 
Tau jet, in the range of 6-15 km s _1 , are in agreement with 
the values predicted by the MCA models, for an assumed 
Mja/Macc ~ 0.1 (which is a bona fide value for these young 
stellar systems and, in particular, for the DG Tau system). This 
result, if combined with the fact that the circunstelar disk that 
surrounds DG Tau shows a velocity pattern [extracted from the 
13 CO(2-l) data by Testi et al. 120021 that is consistent with a 
disk rotating in the same direction as the jet, reinforces the 
MCA scenario. Applying a technique very similar to that of 



1 There are, in fact, some points in their data that do not follow such 
a trend, but this might be due to either small fluxes or the proximity to 
the bow shock of the DG Tau jet; see the discussion in Bacciotti et al. 



Bacciotti et al. |2D5'2l, Woitas et al. (120051 have re-analyzed the 
data from their HST/STIS observations of the RW Aur bipolar 
flow (see Woitas et al. 2002 ), and they found the same trend in 
the radial velocity shifts. 

Coffey et al. ( 120041 . also using the HST/STIS spectrograph, 
found that the bipolar jets from the T Tauri stars TH28 and RW 
Aur (see also Woitas et al. 2005 for the latter source), as well as 
the blue-shifted jet from LkHa321 also show systematic veloc- 
ity asymmetries which are consistent with the results obtained 
by Bacciotti et al. (2002) for DG Tau. In other words, they also 
interpret the radial negative (or positive) velocity shifts as a sig- 
nature of rotation (clockwise or counter-clockwise, depending 
on the sign of the radial velocity shift; see Coffey et al. 2004 
for details). In particular, they found that in the TH28 and RW 
Aur flows, which are bipolar, both the jet and the counter-jet 
appear to rotate in the same direction, and they use this fact as 
an argument in favour of the MCA models. 

These observational data and, more importantly, their inter- 
pretation as a jet rotation signature, lead some authors to use 
them to constrain the wind launching region on the surface of 
the accretion disk. Pesenti et al. (2004) used their self-similar 
MHD disk wind solutions to investigate the possible influence 
of projection or excitation gradients in the relation between the 
rotation signature and the true azimuthal velocity profile. They 
computed cold and warm models with different parameters and 
found that warm models, with an outer radius 2 ~ 3 AU, better 
reproduce the observations of DG Tau. They found that both 
cold and warm disk solutions predicts increasing values for 
the radial velocity shifts for increasing distances from the jet 
axis. However, the values predicted for the velocity shifts in 
the warm solution are systematically smaller than those pre- 
dicted by the cold disk solutions, in better agreement with the 
observations. They also argue that the shift in the radial veloc- 
ity for the computed models are systematically smaller than the 
values given by the true rotational profiles, and that this effect 
is more pronounced near the jet axis. 

In a somewhat different approach, Cerqueira & de Gouveia 
Dal Pino (2004) have performed three-dimensional SPH nu- 
merical simulations of rotating jets. They were able to repro- 
duce the trend in the shift of the radial velocities observed in 
DG Tau by Bacciotti et al. ( 2002 1, with a model of a precessing 
(half-opening angle of 5°), variable velocity (300 + 100 km s -1 , 
sinusoidal velocity variation with a t p „i = 8 years period) jet. 
Cerqueira & de Gouveia Dal Pino (2004) have not computed 
the emission line profiles, that would allow a direct comparison 
between the numerical models and the observed data. These au- 
thors claim that a magnetic field of intensity of the order of « 1 
juG is necessary in order to maintain jet stability against lateral 
expansion due to centrifugal forces. Actually, very recent mod- 
eling of observational data carried out by Woitas et al. (2005) 
for both lobes of the bipolar outflow from RW Aur also suggest 
a magnetic collimation mechanism for the jet. 

These observations and models are increasingly providing 
support for a scenario in which the MCA model matches the 
available observational data (considering the errors estimated 

2 This is the outer radius of the disk region responsible for the 
launching of the wind. 
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from both). However, apart from the very low number of ob- 
served sources with this kind of radial velocity pattern (four 
until now), there are still a few points that should be carefully 
verified before excluding alternative views, which can also ex- 
plain the observed data. For example, Soker ( 2005 ) proposed 
that the interaction between the jet with a warped accretion 
disk can explain the asymmetry in the radial velocity distribu- 
tion across the jet beam, although his model is in need of direct 
observational support (since an inclination between the jet and 
the outer parts of the accretion disk is assumed). 

In this paper, we present the results of a set of fully three- 
dimensional numerical simulations of a non-magnetic 3 jet that 
is allowed to precess and rotate. This, together with an ejec- 
tion velocity time-variability (which produces multiple inter- 
nal working surfaces), represents a complete scenario for the 
Herbig-Haro jets, which are believed to present the observa- 
tional signature of rotation. With a proper calculation of the 
emission coefficients, the cooling of the gas, and the evolution 
of the atomic/ionic species, we are able to produce velocity 
channel maps, and then reconstruct the line profiles for several 
of the observed forbidden and permitted (namely Ha) emission 
lines. We then make a comparison between the radial velocities 
extracted from these computed line profiles and those from pre- 
viously published data of HH jets. 

The paper is organized as follows. In §2, we describe the 
numerical technique and the physical parameters of the simu- 
lations. The results are presented in §3, and in §4 we discuss 
these results and their relation with the available observational 
data. 

2. The numerical method and the simulated 
models 

2.1. The Yguazu-a code 

The 3-D numerical simulations have been carried out with the 
Yguazu-a adaptive grid code (which is described in detail by 
Raga et al. 2000 and Raga et al. 2002 1, using a 5-level binary 
adaptive grid. The Yguazu-a code integrates the gas-dynamic 
equations (employing the "flux vector splitting" scheme of van 
Leer I1982I together with a system of rate equations for the 
atomic/ionic species: HI, HII, Hel, Hell, Helll, CII, CIII, CIV, 
NI, Nil, NIII, OI, Oil, OIII, OIV, SII and SIII. With these rate 
equations, a non-equilibrium cooling function is computed. 
The reaction and cooling rates are described in detail by Raga 
et al. ( 1277021 . 

This code has been extensively employed for simulating 
several astrophysical scenarios such as jets (Masciadri et al. 
120021 Velazquez, Riera & Raga|2004 Raga et al. 120041 Raga 
et al. 120021 1. interacting winds (Gonzalez et al 20(3a| I2004bl 
Raga et al. 2001 a 1 and supernova remnants (Velazquez et 2004 
Velazquez et al. 2001ai. It was also tested with laser generated 
plasma laboratory experiments (Raga et al. 2001b, Sobral et al. 
120001 Velazquez et al. l2001b> . 

The computational domain has dimensions of 1.1 x 10 16 
cm along the x— and y-axes, and of 4.4 x 10 16 cm along the 

3 The discussion of the magnetized case will be left to a forthcom- 
ing paper (De Colle et al. 2005; in preparation). 



Table 1. The simulated models. 



Model 


T pul. Tpre. 


9 


rotation? 


Ml 


8 - 




no 


M2 


8 8 


5 


no 


M3 


8 - 




yes 


M4 


8 8 


5 


yes 



z-axis. The grid has a maximum resolution of 8.59 x 10 13 
cm along the three axes. Assuming a distance of 140 pc (e.g., 
Kenyon, Dobrzycka, & Hartmann 19941, this maximum reso- 
lution corresponds to an angular size of * 0V04. 

2.2. The models 

The parameters were chosen in order to model the observa- 
tions of the DG Tau jet. All the simulated models have a time- 
dependent ejection velocity given by (Raga et al. l2001ct : 



Vo 



, 2tt 
1 + A sin I 1 

Tpul 



(1) 



where vo = 300 km s , and A = 0.33 (giving a jet veloc- 
ity variability in the range of 200-400 km s~') and T pu i = 8 
years. These parameters are appropriate for producing the ve- 
locity structure of the DG Tau microjet (Bacciotti et al. 2002; 
Raga et al. l2001cl Pyo et al.|2003 1. There is also evidence that 
the DG Tau microjet precesses, with a precessing angle of « 5° 
and a precessing period of « 8 years (e.g.; Lavalley-Fouquet et 
al. 2000; Dougados et al. 2000l. In Table^we give the model 
label (first column), the period (in years) for the velocity vari- 
ability and for the precession of the jet axis (second and third 
column, respectively), the half-opening angle of the precession 
cone (fourth column) and, finally, the presence (or not) of ro- 
tation in a given model, expressed by a toroidal velocity pro- 
file (see below). The physical parameters of the models pre- 
sented in Table ^ are, thus, suitable for the study of the DG 
Tau jet. Finally, we choose an initial number density for the 
jet rij = 1000 crrT 3 (with 77 = rij/n e = 10, where n e is the 
initial number density for the environment), and a temperature 
Tj = 10 4 K (with k = TjlT e = 10). We also assume that the 
initial jet H ionization fraction is of 10%. 

In order to investigate the signatures of rotation on the 
line profiles, we have simulated models with and without an 
initially imposed toroidal velocity profile (see Table The 
adopted toroidal velocity is the same as in Cerqueira & de 
GouveiaDal Pino (2004): 



= 8km s 



R 



(2) 



where v* is the toroidal velocity in km s _1 , Rj is the jet radius 
(5.6 x 10 14 cm in physical units) and R is the cylindrical radius. 
We note that this profile has been truncated at a radial distance 
= 0A5Rj, where the toroidal velocity attains a constant value 
of V0 55 km s . 

We should note that the radial dependence of the toroidal 
velocity depends on the details of the MCA wind model that 
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is considered. For example, the model of Pesenti et al. (2004) 
predicts a oc Rr l l 2 law. In the present paper, we restrict our- 
selves to the toroidal velocity law given by equation 0, and do 
not explore other possibilities. 

From the temperature, and atomic/ionic/electronic number 
densities computed in the numerical time-integration, we cal- 
culate the emission line coefficients of a set of permitted and 
forbidden emission lines. We compute the Balmer Ho- line con- 
sidering the contributions from the recombination cascade and 
from n — 1 — > 3 collisional excitations. The forbidden lines 
are all computed by solving 5-level atom problems, using the 
parameters of Mendoza ( 1983 1. 

3. Results 

Figure ^ shows the temporal evolution for the emission maps 
of models M2 (Fig^, b) and M4 (FigQJ, d), assuming an in- 
clination angle with respect to the line of sight of 45° (which is 
the estimated inclination angle of DG Tau; see Pyo et al. l2003i . 
Model M2 presents a variability in both direction and ejection 
velocity, and model M4 also has an initially imposed toroidal 
velocity pattern (see §2). The computed emission lines are: 
Ha and [S II]/16716 (left panels), [O I]/16300 and [N II]^6548 
(right panels). The small precession angle (see Tabled in both 
models produces a gentle wiggle along the jet beam, an effect 
that can be easily seen in Figure ^ The internal working sur- 
faces are morphologically similar in both models. 

The major difference between the M2 and M4 models lies 
in the morphology of the jet head. In the rotating jet model 
M4, the jet head appears to be much less collimated, an effect 
that is due to the centrifugal force that increases the jet opening 
angle. This effect is clearly seen through a direct comparison 
between the emission maps of models M4 and M2 in the head 
region. For model M2, Figures and b shows a collimated 
jet head, while the jet head of model M4 in Figures and d 
shows a clumpy, fragmented structure. Furthermore, there is a 
little deceleration of the jet head in model M4 with respect to 
model M2, an effect that is associated with the de-collimation 
of the jet head. 

For each evolutionary time, and for all of the computed 
emission lines (namely, Ha, [S II]^6716, [O I]^6300 and [N 
II]/16548), we compute velocity channel maps in the -400 to 
100 km s range (with a resolution of 10 km s~'), assuming 
an inclination angle of 45° towards the observer. These veloc- 
ity channel maps allow us to reconstruct the line profiles, since 
we have the intensity of the line for a given radial velocity as a 
function of position in the plane of sky (see §2). Thus, we use 
this datacube to build line profiles for selected regions in our 
computational domain. In particular, we use 3x3 px 2 region to 
define a single, artificial slit (which, for an assumed distance 
of 140 pc and an assumed inclination angle with respect to the 
line of sight of 45°, gives us a « 071 x 071 slit). An intensity is 
determined for each velocity channel, as a function of the slit 
position in the computational domain. These intensity points, at 
10 km s radial velocity intervals, are then combined to give 
the line profile for each slit (as a function of position). 

We define four regions along the jet axis, namely regions 
I, II, III and IV, at distances from the jet inlet of O'.'l, 072, 073 



and 074, respectively (assuming a distance of 140 pc and an 
inclination angle of 45°; see §2). For each one of these regions, 
we construct seven 07 1 x07 1 slits (as described above), that are 
placed side by side parallel to the jet axis with one of them lying 
exactly on the jet axis (in order to have slits symmetrically dis- 
posed on each side of the jet). This configuration (seven aligned 
slits across the jet beam and in four regions along the jet beam) 
was chosen in order to mimic the observational conditions de- 
scribed in Bacciotti et al. (2002i. In Figure[2]we show the slits 
and their positions on the plane of the sky. We note that, in fact, 
there is an overlap of 1 pixel {k 0707) between the slits in each 
region (see Figure^)- 

In Figure^, we show the Ha line profiles for region IV of 
model M4 at t — 10 yr. We note the presence of at least three 
components in several of the line profiles shown in Figure |3| 
The emission in the radial velocity range from » to « -100 
km s _1 is related to the leading bow shock, which surrounds the 
jet beam. We assume that the emission in the interval from -300 
km s _I to -100 km s _1 is related to two components, namely, 
the high- and medium-velocity components (hereafter, HVC 
and MVC, respectively). We have performed standard mini- 
mum ^-Gaussian fits to obtain the MVC, for the seven slits in 
each of the four regions defined in Figure [2] Figure |5j3 shows 
the data (full line) combined with the gaussian fit (crosses) for 
the whole radial velocity interval (left) and for the radial veloc- 
ity interval ranging from -300 km s _1 to -100 km s _1 (right). We 
also show the gaussian fits decomposing the HVC and MVC. 
We have performed these fits for all models in Tabled and for 
different evolutionary times. This procedure allows us to inves- 
tigate both the behaviour of the MVC for each emission line as 
well as its dependence with the imposed initial conditions. 

Two distinct evolutionary times were chosen in order to 
investigate the influence of the presence of an internal work- 
ing surface (IWS) close to the region where the the slits were 
placed. In particular, at t — 10 years (see Fig.Q, an IWS is 
just being formed (close to the jet inlet, in regions I-IV) due 
to the chosen T pu i - 8 years. On the other hand, for t = 15 
years, the regions I-IV are spatially detached from the position 
of the nearest IWS (see Fig.|2}. In Figure|4]we show the radial 
velocity of the MVC for regions I, II, III and IV (from top to 
bottom) as a function of distance from the jet axis (from -072 
to 072; the in the abscissa corresponding to the position of 
the jet axis) for t = 10 years (left) and t — 15 years (right). 
The radial velocity of the MVC was extracted from the profiles 
of the Ha, [O I]^6300, [S II]/16716 and [N II] ,16548 emission 
lines, using the procedure described above. Although the MVC 
from different emission lines shows slightly different behav- 
ior, the distribution of the radial velocity on both sides of the 
jet axis for a given emission line is perfectly symmetric for the 
Ml model, regardless the analyzed region or evolutionary time, 
as we can see clearly in Figure^] 

In Figure |5] we show the radial velocity (V ra d', left) for 
regions I, II, III and IV (from top to bottom) of model M2 
(with precession but with no rotation) as a function of dis- 
tance from the jet axis (from -073 to 073; as in the previ- 
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ous diagram ), as well as the difference (AV ra d) between the 
radial velocity taken at symmetrical positions with respect to 
the jet axis (right), at t — 10 years (top diagrams) and t = 15 
years (bottom diagrams). The precession of the jet axis intro- 
duces an asymmetry in the radial velocities as a function of 
distance from the jet axis, as we can see in the left diagrams 
of Figure 13 This effect is more pronounced at t = 10 years, 
as indicated by the difference in the radial velocities, AV ra d, in 
Fig. |5] (top right diagram), and this is due to the presence of 
an IWS close to regions I-IV. However, such an asymmetry in 
the distribution of the radial velocity is still present at t — 15 
years. As we can see in Figure |5] (bottom diagrams), AV ra d 
is systematically negative for most of the computed emission 
lines. Although we cannot establish an unique behaviour for 
all of the computed lines (considering the distinct regions and 
evolutionary times), we note that some of them show a clear 
trend of smaller AV ra d for smaller distances from the jet axis 

(V ra d,3 ~ V ra d,5 < Vrad,2 ~ V ra d,6 < Vrad,\ ~ Vrad.l)- 

In Figure |6] we again show the radial velocity (V ra ,/; left) 
as a function of distance from the jet axis for regions I, II, III 
and IV (from top to bottom), as well as the difference (AV r „d) 
between the radial velocity taken at symmetrical positions with 
respect to the jet axis (right), at t — 10 years (top diagrams) and 
t — 15 years (bottom diagrams), but now for model M3. This 
model presents variability in the velocity of injection as well 
as an initially imposed rotational profile inside the jet beam, 
but has no precession (see §2 and Table The asymmetry 
in the radial velocity distribution, with respect to the jet axis, 
can be clearly seen in Fig. [6] In particular, AV m d is negative 
for almost all emission lines, regions and evolutionary stages 
(with the exception of the [N II]/16548 emission line, at t — 10 
years in region II). However, in opposition to model M2, AV ra d 
appears to be greater for smaller distances from the jet axis (at 
least for some emission lines). 

We present the radial velocity distribution for model M4 
in Figure Model M4 has variability in both direction and 
velocity of injection, as well as an initially imposed rotation 
profile inside the jet beam (see Table [Q. The asymmetry in 
the radial velocity distribution on both sides of the jet axis is 
more pronounced than in models M2 and M3. This is clearly 
seen from the radial velocity difference diagrams (right panels 
in Figure 0, which show AV n ,d of the order of up to « -80 
km s _1 . As in model M2, there is a general trend in the AV ra d 
data that shows that it increases as a function of distance from 
the jet axis, although such a trend is not followed by the He, [O 
I]/16300 and [S II]/16716 emission lines in region IV at t — 10 
years and region I at t — 15 years (see the right diagrams of 
Figure 0. 

In Figure [8] we show the radial velocity shift as a func- 
tion of distance from the jet axis for regions I to IV (from 

4 We note that the precession introduces a small shift on the jet axis 
regarding the central pixel in the x-direction in our simulations, i.e., 
x=65. However, due to the very small precession angle chosen (see 
Table 0, the displacement of the slits with respect to the jet axis is 
of the order of < (X'02. This small value should introduce negligible 
errors in the radial velocity determinations. Moreover, this value is 
quantitatively comparable with the deviation from the jet axis in ob- 
servations of Bacciotti et al; see Figure 8 in Bacciotti et al. 1 2002 1. 



top to bottom), and for models Ml (crosses), M2 (asterisks), 
M3 (diamonds) and M4 (squares), for t — 10 years (left) and 
t — 15 years (right). Each point in this diagram corresponds 
to a mean value for the radial velocity shifts, assuming that 
AV m d — Z/Li AV ra dj/4; where the different values of / indi- 
cate the Ha, [O Ij/16300, [N II]i6548 and [S 11^6716 emis- 
sion lines. We clearly see a general trend of increasing radial 
velocity shifts for larger distances from the jet axis, followed 
by all the models (with the expected exception of model Ml, 
which shows a null AV, as discussed previously). Such a trend 
is similar to the one discussed in Pesenti et al. (2004i for their 
disk solution model (and also followed in part by the data of 
Bacciotti et al l2002> . 

4. Discussion and conclusions 

We have presented a set of 3D numerical simulations of vari- 
able velocity jets, from which we have obtained detailed pre- 
dictions of the emission line profiles of a set of emission lines 
(Ha, [O I]/16300, [S II]^6716 and [N II]^6548). The models in- 
clude an axisymmetric jet, a precessing jet, a rotating jet, and, 
finally a jet with both precession and rotation (see Table 1 and 
section 2). 

Following Bacciotti et al. (2002), we carry out double- 
Gaussian fits to the line profiles, and then study the behaviour 
of the radial velocity of the MVC (medium velocity compo- 
nent) as a function of distance from the jet axis, for several 
positions along the jet flow. We find that while the axisymmet- 
ric model of course predicts radial velocities whith a symmetric 
behaviour on both sides of the jet axis, the other three models 
predict asymmetric behaviours. 

The models with precession and with rotation predict sys- 
tematic offsets of |AV ra d| < 20 km s -1 for the MVC between 
symmetric positions on both sides of the jet axis. These sys- 
tematic offsets are seen at most positions along the jet axis, and 
for many of the emission lines. The magnitudes of these off- 
sets are similar to the ones observed (for many positions and 
emission lines) by Bacciotti et al. (2002), Coffey et al. (2004) 
and Woitas et al. ( 2005 ). In the model with both precession and 
rotation (in which the precession was assumed to have a pro- 
grade direction), the two effects combine to give larger offsets 
of- 40- 80 km s _1 . 

We have also computed a model with rotation and a retro- 
grade precession. This model (not presented in this paper) pro- 
duces side-to-side offsets in the MVC radial velocity which are 
much smaller, and which do not show clear trends as a function 
of distance from the jet axis. 

From these results, we conclude that the radial velocity 
asymmetries observed in the DG Tau, TH28, RW Aur and 
LkHa321 microjets (see Bacciotti et al. 2002; Coffey et al. 
2004; Woitas et al. I2005i can be modeled either as due to a 
rotation of the gas in the initial jet cross section, or as due to 
a precession of the ejection direction. Interestingly, though a 
precession results in a point-symmetry in the "precession spi- 
ral" morphology of the jet beam, the direction of the resulting 
rotation of the ejected material is of course the same one for 
both the jet and the counterjet. Also, a "precession spiral" mor- 
phology can be produced by an orbital motion of the source 
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around a binary companion (Masciadri & Raga 2002 Gonzalez 
& Raga 2004 1, resulting in similar radial velocity patterns, but 
with a morphology with mirror symmetry between the jet and 
the counterjet. 

It is therefore unclear at this time whether the systematic ra- 
dial velocity asymmetries observed by Bacciotti et al. (|2002 1, 
Coffey et al. ( 2004 1 and Woitas et al. ( 2005 ) in the cross sec- 
tions of several microjets are due to jet rotation or precession 
(or to a combination of both effects). In order to proceed fur- 
ther, it will be necessary to proceed with more detailed mod- 
elling, in which not only the radial velocity asymmetries are 
considered. A study which also considers the locci of the jet 
and counterjet on the plane of the sky, resulting in models 
which simultaneously reproduce the morphology and the radial 
velocity structure of the observed outflows, might be able to 
solve this problem. Of the four microjets observed by Bacciotti 
et al. d200lt . Coffey et al. ( 120041 and Woitas et al. (2005 1, the 
jet from RW Aurigae shows the chain of knots with the best 
alignment. Because of this lack of evidence for a precession 
spiral structure in the jet beam, the RW Aur jet is probably the 
most certain case for an observed rotation of the jet beam. 

Interestingly, the referee has pointed out to us the very re- 
cent results of Cabrit et al. (2005 1, who detect a disk around 
RW Aur rotating in the opposite sense of the "jet rotation" seen 
by Woitas et al. (2005 1. This result in principle eliminates a jet 
rotation interpretation for the side-to-side radial velocity asym- 
metries observed in the RW Aur jet. Therefore, the remaining 
possibilities are that the jet is precessing (producing the asym- 
metries explored above), or that it has another intrinsic side-to- 
side asymmetry (see, e. g., Soker 2005 and Cabrit et al. 2005). 
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Fig. 1. a) Temporal evolution (f = 10, 15, 20 and 25 years, from left to right) of the emission maps for the Ha (top) and [S 
II]i6716 (bottom) lines, computed for model M2 (see Table[TJ. b) The same as in a) for the [O I]^6300 (top) and [N II]/16548 
(bottom) lines, c) Temporal evolution (f = 10, 15, 20 and 25 years, from left to right) of the emission maps for the Ha (top) and 
[S II]^6716 (bottom) lines, computed for model M4 (see TableQ}- d) The same as in c) for the [O I]^6300 (top) and [N II]/16548 
(bottom) lines. The jet in each map was assumed to have an angle with respect to the line of sight of 45°. The bars on top of each 
diagram give the normalized flux. [See the electronic edition of Astronomy & Astrophysics for a color version of this figure.] 
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Fig. 2. left) An [O II]/16300 line emission map for the M4 model in a region near to the jet inlet, at t=15 years. The coordinates 
are in units of pixels, and the scale bar indicates the normalized intensity, right) A closer view of the jet base, corresponding 
to the region indicated by the rectangular box in the left-panel, indicating also the artificial slit positions. Regions I, II, III and 
IV are successively more distant from the jet inlet, while positions 1, 2, 3, 4, 5, 6 and 7 corresponds to different position across 
the jet radius. Note that each slit have a superposition of 1 px in the direction transverse to the jet axis. For the sake of clarity, 
we show in this figure only a few (one or two) of the seven slits for each region. See the text for discussion. [See the electronic 
edition of Astronomy & Astrophysics for a color version of this figure.] 



10 



Cerqueira et al: Emission lines from rotating jets 



3.0 
2.5 
2.0 
1.5 
1.0 
0.5 

3.0 
2.5 
2.0 
1.5 
1.0 
0.5 

3.0 
2.5 
2.0 
1.5 
1.0 
0.5 



slit 1 




slit 2 




slit 3- 




slit 4 



slit 5 




slit 6: 




7 



-300-200-100 



(km s ') 



300-200- 
V rod (kr 

M4 



100 
i S-) 




-400-300-200-100 



100 



Mode 
t= 1 years 
Ha 

Region IV 



(km s"') 




Fig. 3. a) Ha line profile for model M4, at t — 10 years, for region 4 (see Figure|5]for slit and region definitions). The slit number 
is given on the top-right corner of each diagram. The intensity is in units of 10~ 6 erg cirT 2 s~' . b) Left: The Ha line profile of slit 
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Fig. 4. Radial velocity as a function of distance from the jet axis for regions I to IV (from top to bottom), for t — 10 years (left) 
and t — 15 years (right) and model Ml (see Table^. The distance is in units of arcsec, considering a distance of 140 pc from the 
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. 6. The same as in Figure|5J but for the M3 model (see Tableful. 
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. 7. The same as in Figure|5J but for the M4 model (see Tableful. 
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Fig. 8. The radial shift as a function of the distance from the jet axis for regions I, II, III and IV (from top to bottom), for models 
Ml (crosses), M2 (asterisks), M3 (diamonds) and M4 (squares), for t = 10 years (left) and t = 15 years (right). 



